The WIMP capture process for dark stars in the early universe 
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The first stars to form in tiie universe may iiave been dark stars, powered by dark matter anniliilation 
instead of nuclear fusion. The initial amount of dark matter gathered by the star gravitationally can 
sustain it only for a limited period of time. It has been suggested that capture of additional dark 
matter from the environment can prolong the dark star phase even to the present day. Here we show 
that this capture process is ineffective to prolong the life of the first generation of dark stars. We 
construct a Monte-Carlo simulation that follows each Weakly Interacting Massive Particle (WIMP) 
in the dark matter halo as its orbit responds to the formation and evolution of the dark star, as it 
scatters off the star's nuclei, and as it annihilates inside the star. A rapid depletion of the WIMPs 
on orbits that cross the star causes the demise of the first generation of dark stars. We suggest 
that a second generation of dark stars may in principle survive much longer through capture. We 
comment on the effect of relaxing our assumptions. 



The first stars to form in the universe (Population III 
stars) provide the light that ends the dark ages, con- 
tribute to reionization, inject metals into the interstel- 
lar medium for later stellar generations, and might give 
the seeds for the massive early black holes observed in 
quasars. All these effects depend on the properties of the 
first stars, which in turn depend strongly on their masses. 

Dark matter can play a crucial role in the formation 
of the first stars and in particular can change their 
masses and astrophysical implications dramatically 
The first stars form at the centers of early dark mat- 
ter halos, wlien the universe was young and much more 
dense than today. Not only were the first stars formed 
in regions of very high dark matter densities, but they 
also dragged in more dark matter gravitationally while 
they were forming (a process often called "adiabatic con- 
traction" although it may not occur adiabatically) . If 
dark matter consists of Weakly Interacting Massive Par- 
ticles (WIMPs) that annihilate among themselves, their 
annihilation could provide enough energy to stop or de- 
lay the formation of Pop III stars, giving rise to a new 
stellar phase called a Dark Star. This annihilation is the 
same physical process that creates WIMPs thermally in 
the early universe and naturally gives the correct cos- 
mic density for WIMPs to account for the observed dark 
matter. 

Once the initial amount of dark matter accumulated 
via adiabatic contraction runs out, the dark star phase 
ends and the star contracts onto the zero-age zero- 
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metallicity main sequence [2]. It has been suggested 0-01 
that capture of additional dark matter from the environ- 
ment around the star can prolong the dark star phase 
even to the present day. 

Here we show that this capture process is ineffective to 
prolong the life of the first generation of dark stars. For 
these stars, the great majority of the dark matter around 
them is accumulated by adiabatic contraction at the time 
of formation and ends up in orbits that are bound to 
the dark star. Part of this dark matter is captured by 
the dark star once the gas density in the star is high 
enough. Essentially all of the dark matter in the dark 
star annihilates away during its contraction towards the 
main sequence at the end of the dark star phase. The 
rest of the dark matter remains close but outside the 
star and orbits around it without crossing it. We call 
this dark matter population the dark matter envelope. 
Once the dark matter crossing the star has been depleted, 
there is no additional dark matter population that can 
be efficiently captured. This ends the first generation of 
dark stars. 

Notice that in this paper we have extended the mean- 
ing of the word capture in the following sense. By cap- 
ture, we mean the concept of bringing WIMPs from "out- 
side" to "inside" the star through energy loss in WIMP- 
nucleon scattering, which increases the fraction of time 
a WIMP spends inside the star. Originally capture was 
defined as the scattering of dark matter particles in the 
Galactic population off individual nucleons in the Sun 
and into bound orbits 0, 01 • In the case of Pop III dark 
stars, most of the WIMPs accumulated around the star 
by adiabatic contraction are bound to the star, and there 
is no relevant unbound population. Hence the original 
formalism for capture, as used in 0-13 1 does not apply. 
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We instead define capture as the scattering of dark mat- 
ter particles off nucleons in the star without regard to the 
bound or unbound nature of their orbits. This definition 
allows for a clear separation of scattering and annihila- 
tion processes. 



I. INTRODUCTION 

WIMPs constitute a negligible fraction of the mass of 
the dark star (less than a few percent), but the WIMP 
density is high enough for WIMP annihilation to be- 
come the dominant power source. This was proposed 
in Ref. [1], which introduced the term Dark Stars to re- 
fer to Pop III stars in a stellar evolution phase powered 
by dark matter. Later the term was extended in Ref. ^ 
to any star powered by dark matter annihilation, previ- 
ously called a dark matter burne r I9l4l3l . Work on Pop 
III dark stars include Refs. [ll-H, 1 13142 ill : work on other 
kinds of dark stars include Refs. |8l-[l3. [22 - [25| . In this 
paper, we consider Pop III dark stars exclusively. 

In dark stars, the dark matter becomes the dominant 
heat source despite its comparably low density because 
WIMP annihilation converts mass into energy much more 
efficiently than nuclear fusion. As in Ref. we assume 
that 2/3 of the WIMP annihilation energy is deposited as 
heat in the star and the other 1/3 is lost, predominantly 
to neutrinos. 

As the dark star forms, the deepening of the gravita- 
tional potential boosts the halo's central WIMP density 
in the proximity of the star (adiabatic contraction) . An- 
nihilation within the WIMP population heats up the dark 
star, but also decreases the dark matter density. How- 
ever, as long as the dark star keeps gaining mass and 
increases its density, the central dark matter population 
inside the star is replenished by the increased gravita- 
tional pull. The energy injected by WIMP annihilation 
halts the contraction of the star, keeping the outer re- 
gions cool for a time longer than a normal Pop III star 
without dark matter. The dark star is then supported by 
the energy from the annihilation of WIMPs. Initially, it 
is a low-density giant star with low surface temperature 
[l|. Later, by accreting mass from the surrounding gas, it 
grows to a mass of several hundreds of solar masses, still 
remaining a low-density giant star (radius of thousands 
of solar radii). The dark matter has given the star more 
time to accrete mass, making the first stars more massive 
than without the help from dark matter annihilation. 

As the dark star becomes bigger and more massive, 
more energy is needed to support it. To sustain the star, 
dark matter needs to be supplied at a higher rate. Even- 
tually the adiabatic contraction can no longer keep up 
with the dark star's increasing demand of dark matter. 
At this point, the dark star starts to contract since the 
dark matter supply rate is no longer enough to support 
it. When the volume of the dark star decreases, the total 
amount of energy supplied by dark matter annihilation 
decreases with the volume. After the contraction, the 



star approaches the main sequence. When nuclear fu- 
sion starts, the radiation feedback halts the accretion of 
baryons, putting an end to the star's growth in mass. 
Then the dark matter can no longer be gravitationally 
replenished by adiabatic contraction, putting an end to 
the initial dark star phase. The end product is a normal, 
although very massive, main sequence star in the sense 
that it produces energy through nuclear fusion and no 
longer through WIMP annihilation. 

It was suggested in Refs. that the WIMP popula- 
tion at the end of the dark star phase may be replenished 
by capture of WIMPs from the surrounding dark matter 
halo. WIMPs passing through the star can loose energy 
in scattering off nucleons in the star. If a WIMP scatters 
many times it will eventually loose enough energy to sink 
to the center and contribute to the WIMP density inside 
the star. If enough energy is injected through the annihi- 
lation of the captured dark matter, the nuclear burning 
in the star is slowed down, possibly to essentially a full 
stop. This object would then be powered by dark matter 
and, as long as the dark matter fuel is constantly refilled, 
the lifetime of the star is substantially prolonged, possi- 
bly to time scales comparable to the age of the universe 
("eternal dark stars" 041). 

For WIMP capture to prolong the life of a dark star, it 
is essential that the dark matter is efficiently captured via 
WIMP-nucleon scattering. Refs. estimated the rate 
of WIMP capture using a formula by Gould 0|- How- 
ever, Gould's formula assumes a constant flow of WIMPs 
coming onto the star from far away. It also assumes that 
the WIMP population is Maxwellian and is not bound 
to the star. However, a Pop III dark star forms thanks 
to the adiabatic contraction of its own dark matter halo, 
and the WIMPs responsible for the high concentration 
of dark matter in and around the star are gravitation- 
ally bound to the star. Thus also the WIMPs accessible 
for capture are gravitationally bound to the star, and 
Gould's formula does not apply. 

In this work, we investigate in detail the star's abil- 
ity to access dark matter from the halo it formed in. 
For this purpose, we construct a Monte-Carlo simulation 
that follows the halo WIMPs as their orbits respond to 
the formation of the star, as they scatter off the star's 
nuclei, and as they annihilate inside the star. Since the 
WIMPs are fully traced, not only can we analyze their 
annihilation and scattering behavior, but also their phase 
space distribution around the formed star. 

We find that WIMP capture is not important until the 
dark star contracts on its way to the main sequence and 
becomes dense enough to be a viable dark matter target. 
We also find that WIMP capture lasts a relatively short 
amount of time as the WIMPs available for capture are 
quickly depleted and their orbits are not replenished. 

Throughout this paper we use as fiducial values a 
WIMP mass — 100 GeV and an annihilation cross 
section times relative velocity aannV = 3 x 10~^^ cm'^/s. 
For scattering off nuclei, since the spin dependent scatter- 
ing cross section is far less constrained than the spin inde- 
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pendent cross section, we look at the case of a large spin- 
dependent scattering cross section CTgcatt — 10^'^^ cm^, 
as well as at a zero scattering cross section (absence of 
scattering). 

Section II describes our Monte-Carlo (initial phase 
space distribution, adiabatic contraction, capture and an- 
nihilation) . There we also present a new quadrature for- 
mula for the phase-space distribution of a Navarro-Frenk- 
White halo. Section III presents our results, which 
are then discussed in Section IV. 



II. METHOD 

The motivation for this work is to find out how the 
dark matter in the NFW halo responds to the change in 
gravitational potential as the star forms, and what then 
happens to the dark matter with time. Following the 
choices of Ref. [l| we investigate the case of a single star 
forming at redshift z — 20 in the center of a primordial 
halo of mass IO^Mq with a Navarro, Frenk and White 
(NFW) profile [H consisting of 15% baryons and 85% 
dark matter. 

The Monte-Carlo simulation picks individual WIMPs 
from the original NFW halo and follows these WIMPs 
and how they respond to the changing gravitational po- 
tential as the star forms. Through picking and studying 
WIMPs in phase space the full orbit of the WIMP inves- 
tigated is known. With this information the adiabatic 
invariants of the orbit can be determined and with that, 
as long as the star formation process is slow enough to 
be adiabatic, the new orbits of the WIMPs can be found 
through all stages of the stellar evolution. The phase 
space distribution of the dark matter as the star evolves 
is thus followed. 

WIMPs spending time inside the dark star can anni- 
hilate with other WIMPs or they can scatter off nuclei 
in the star. If the WIMPs annihilate, their energy (mi- 
nus the energy lost to neutrinos) is released to heat the 
gas. If the WIMPs scatter before having time to anni- 
hilate, the energy loss and then the new WIMP orbit 
is calculated. The new orbit of the scattered WIMP is 
then followed and the WIMP is again given the chance 
to scatter or annihilate. As the star evolves, the future 
orbits of the WIMPs which have not annihilated are cal- 
culated, so that the WIMPs are given new chances to 
scatter and/or annihilate at the later stages of the stellar 
evolution. Simulating many WIMPs gives the injected 
heat from dark matter into the dark star as a function of 
time, including the effect of WIMP-nucleon scattering. 

The evolution used for the dark star is taken from 
Ref. More precisely, it is their canonical case of a 
100 GeV WIMP with Tan-McKee mass accretion ^ 
and without WIMP capture (i.e. without WIMP-nucleon 
scattering) . To follow the WIMPs through the formation 
of the star, the stellar evolution is divided into 116 steps 
(taken from f^) and the WIMP orbits are found for each 
of these stages. Our work only considers how the dark 



matter halo evolves given the stellar evolution. How the 
dark star responds to the change of the dark matter heat- 
ing from WIMP capture will be taken into account in an 
upcoming publication (28| . 



A. Initial dark matter distribution 

The first step of the Monte-Carlo is to pick WIMPs 
from the original NFW halo, i.e. to find the phase space 
distribution of the original halo. Spherical symmetry is 
assumed to hold throughout, for the star forming in the 
center of the NFW halo and at later times. 

The density profile of an NFW halo is given by 



Sc 



Pc {r/rs){l + r/rsy 



(1) 



where — r2oa/c is a characteristic radius and pc = 
3H^/8ttG is the critical density. Sc is a dimensionless 
concentration parameter given by the requirement that 
the mean density within r2oo is given by 200pc, as in 
Ref. [26|. We take the total halo mass to be 10^ Mq and 
the concentration parameter c = 2. Then 6c — 1235, — 
2.28 X 10^0 cm = 74.5 pc, and pc = 2.5 x lO^^eg/^^^s ^ 
0.014 GeV/cm3. 

For a spherically symmetric density distribution and an 
everywhere isotropic velocity dispersion, the phase-space 
distribution function can be found using Eddington's for- 
mula [29|]: 



1 



d* d^p 1 / dp 



(2) 

Here f{£) has units of mass per unit phase space volume 
(g cm~'^ (cm/s)""^), £ = —E/m^ is the negative of the 
total energy E per unit WIMP mass to^, and \['(r) — 
—4){r) is the negative of the gravitational potential 0(r). 
The condition that the WIMPs are gravitationally bound 
to the NFW halo imphes < f < *(0). For an NFW 
halo, ^'(r) is given by: 



NFW 



(r) = G 



M{r'] 



dr' = 47rG(5cPc 



jln(r/rs 



r/r 



(3) 

We have written M{r) for the halo mass enclosed within 
radius r. 

After a very lengthy but straightforward calculation 
starting from Eddington's formula, the distribution func- 
tion for an NFW profile can be written as: 



/nfw('^^nfw) 



9{r)- 



1 



where 



9{t) 



\/ ^NFW — 'I'nFw('") 

r(rs — 3r -I- 6^) 



idr, 



(4) 
(5) 
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FIG. 1: The phase-space distribution function 
/npw(£nfw) (red solid line) and the fNFW-distribution 
/NPw(fNFw)\/l - fNFw/'I'NFw(r) at r = 10~^rs (blue 
dotted hne) for our 10*' A/q, c = 2 NFW halo. 

with 

e = (r, +r)ln(^l + ^^ . (6) 

Moreover, the lower limit of integration in Eq. (j4|), r^^^, 
is the solution of 5'nfw(''£^^) — ^^nfw- This solution 
cannot be written in compact form. The NFW phase- 
space distribution /nfw(^^nfw) is plotted in Fig.[TJ 

Knowing /nfw('^-nfw)! the Monte-Carlo can pick indi- 
vidual WIMPs from the original halo. In practice, we first 
pick a spatial location r from the distribution 47rr^/9(r)dr, 
then we select a W IMP energy from t he distribu- 
tion 47r/NFw('^^NFw)\/2[^'NFw('') " ^MFwld^NEW, and 
finally we generate a direction for the WIMP velocity 
from a uniform distribution on the sphere. The latter 
two distributions use the isotropy in velocity space. The 
WIMP angular momentum per unit mass J is easily de- 
termined from ^NFW, ^NFw('') and the velocity direction 
at the point r: J = r sin 0-^/2 [5'nfw('") — ^nfw], where 
is the angle between the velocity vector and the outward 
radial direction. With this information, the WIMP's or- 
bit is completely determined up to the spherical symme- 
try. 

There is a solvable numerical difhculty in generating 
^NFW) connected to the d ivergence of J'nfw(^nfw) as 
i^NFW ~^ ^NFw(O). The divergence is in itself not a 
problem because a finite radius r of the WIMP location 
avoids it as £nfw ^ ^'nfw('') < ^'nfw(O); thus in prac- 
tice the divergence point is never reached since M{r) -> 
as r — T' 0. The finite value of r guarantees that the distri- 
bution from whi ch the WIMPs are g enerated, namely 
47r/NFW NFW ) V2[*nfw(?') - ^NFwldfNFW , remains 
finite everywhere. This distribution is shown in Fig. m 
(blue dotted line), rescaled to fit f{£) for small £. 
Even though the distribution from which the WIMPs are 
picked is finite, the divergence in f{£) still produces a 



sharp peak near iPnfw = ^'nfw(''), as seen in Fig. [T] 
In the Monte-Carlo, the peak is typically even sharper 
than in Fig. [1] as r is typically smaller for the interest- 
ing WIMPs, making the distribution numerically chal- 
lenging to use in the Monte-Carlo. The latter issue can 
be treated by going back to Eq. ([4]); Taylor expanding 
around the critical point £nfw — 'I'nfw('')i and inte- 
grating in r, one finds that the peak can be tamed by 
rewriting the velocity distribution in the new variable 

Z = 1/ (*nFw(0) - ^NFW)- 

B. Adiabatic contraction 

As the star forms, the gravitational potential changes 
from «'nfw('^) to «'*(r)-h4'haio(r) + 4'gas(r-), where 4'^(r) 
is the potential generated by the star, \I'haio('") is the po- 
tential generated by the dark matter halo while it adia- 
batically contracts, and '^gas{r) is the potential generated 
by the gas surrounding the star. 

The star influences the potential out to a radius of 
influence rinfl, which can be estimated by equating the 
mass of the star M^, to the amount of dark matter in the 
halo contained within nnfl, i.e. setting Afhaio(nnfi) = M*. 
To find Afhaio(?'infi) one can use adiabatic contraction in 
the circular orbit approximation, in which WIMP or- 
bits do not cross each other so that -MNFw(^hifl^) ~ 
Mhaio(nnfl)- Here r^f^ is the radius the WIMP orbit 
had before it adiabatically contracted to rinfl. Using this 
relation with the NFW density profile we obtain r?^,^^ ~ 
^Ah/{2T:pc5crs) ^ 0.025rs V'IVSOOMq. For circular 
orbits, rAf (r) is an adiabatic invariant. After the forma- 
tion of the dark star, the mass doubles inside a WIMP 
orbit originating at r^^^ and the invariance of rM{r) 
then tells us that ri„fl = rfj^^ /2 ~ O.Olr^ y^M^/SOOM©. 
Another way to estimate the radius of influence is to es- 
timate how far a WIMP would go if it would leave the 
surface of the star at a speed equal to the star's escape ve- 
locity Wesc = \/2GM7JrI. WIMPs which travel further 
out in the halo are not gravitationally bound to the star. 
For the density distribution of the original NFW halo, 
one finds in this way that rf^^ — ^/M^JJ2np^|6^^rs) , 
which is the same expression as in the previous estimate 
of the radius of influence. In conclusion, a good estimate 
of the radius of influence is rjnfl ~ O.Olrs for a final star 
of mass Mi, — 800Mq and radius i?^ = 4 x 10^^ cm. 

Within the star and its immediate surroundings, the 
star contribution ^'*(r) dominates over ^haio(''), and 
thus in that region we neglect ^'haio('')- This approxi- 
mation is good out to the radius of influence of the star 
''infi- By neglecting the halo component ^haio(''), we have 
implicitly assumed that all WIMPs close to the star are 
bounded to the star. In Section IV we show that the 
contribution of unbound WIMPs is negligible. 

Similarly, the density of the gas surrounding the star 
is so low that the potential ^'gas(?') it generates can 
be neglected. For example, for a mass accretion rate 
M = 10~^ Mo/yr, in the middle of the range of the Tan- 
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McKee [27] and O' Shea-Norman [30'| accretion rates, the 
hydrodynamic equation gives a gas density that falls as 
j,-3/2 fj-Qjjj ^jjg surface of the star, where it is approxi- 
mately pga.(i?*) = M/{AttRWGM,R,) - 10-10 g/cm3. 
This is to be compared with an average stellar den- 
sity of ~ Gg/cm'^ for the final stages of the dark star 
(A/* = 8OOAf0, i?* = 4 X 10^1 cm). As a consequence, 
the gas mass Mgas(?') inside a radius r is roughly one half 
of the stellar mass M^, for r = rinti- Hence also 'i>gas{r) 
can be neglected for the region r < rind in which we are 
interested. 

Within the star and its immediate surroundings, the 
star contribution vE'*(r) dominates over both vE'haio('') and 
^'gas('''), and thus in that region we neglect vE'haio(?') and 
^gas(?')- This approximation is good out to the radius of 
influence of the star r-mfi. By neglecting the halo compo- 
nents ^^haioif) and 5'gas(?'), we have implicitly assumed 
that all WIMPs close to the star are bound to the star. 
In Section IV we show that the contribution of unbound 
WIMPs is negligible. This tells us that we can safely as- 
sume that the potential in the region of interest is dom- 
inated by the potential of the star ^'*(r). The potential 
inside the star can be obtained numerically from the tab- 
ulated dark star evolution models in Ref. ^]. Outside the 
star, the potential vE'*(r) is simply Keplerian. 

The individual WIMP's response to the adiabatic con- 
traction of the forming star can be followed using the 
adiabatic invariants 



their orbits are Keplerian. The solution is |3l| 



£ ^ 



{Jr + 271 jy 



(10) 



For WIMPs that cross the star the orbit is not Keple- 
rian and £ needs to be found numerically as the potential 
inside the star is only known numerically. This is done 
by a binary search in f , as J7r is a decreasing function 
of £ at fixed J^, using the circular solution for the given 
J' as an upper limit for £. However, also the integration 
limits Tniin and Tmax depend on £ and have to be com- 
puted with every guess of £. For the orbits of interest 
here, r^im and sometimes also Tmax, are less than i?^ and 
need to be found numerically. This is also done by a bi- 
nary search since the expression inside the square root in 



Eq. ^ is positive for 



and negative for 



T > Tmax and r < Tmin. The starting point for this search 
is the radius Tcirc of the circular orbit for the given JZ, as 



one can show that r. 



mm ^ ^circ ^ ^ max 

2 _ ^3 

circ 



The circular ra- 



dius Tcirc can be found using = -r^;j,^[d^'(rcirc)/dr], 
which can be solved efficiently in the Monte-Carlo since 
the right hand side depends only on the stellar configura- 
tion and does not need to be calculated for every WIMP 
analyzed. To find the orbits of the interesting WIMPs, £ 
needs to be determined for all the different stellar stages 
for which the WIMP orbit crosses the star. 



J and Jr = 2 



y/2[^'{r)-£]-jyr^ dr, (7) 



where is the angular momentum per unit WIMP mass 
and J"r is the radial action per unit WIMP mass. These 
invariants assume the same value before, during, and af- 
ter the formation of the star, as long as the WIMPs do 
not scatter. This allows us to generate the distribution 
of WIMPs as they enter the star for the first time. 

The energy £ is not an adiabatic invariant, however, 
and needs to be recomputed as the potential changes. For 
this, we first compute Jr for the NFW potential using the 
energy £nfw generated from /nfw(-£'NFw) as explained 
above: 



V2[*NFw(r) - fNFw] - jyr^ dr. 



(8) 

As f^j- is an adiabatic invariant J7r = J^^^ ■ We then 
solve the following equation for the new energy £: 



2 / ^2[^4r)-£]-jyr^ dr = X 



NFW 



(9) 



Here the limits of integration r^mi and ' 



are the zeros 



of the square root. In writing this equation, we have 
assumed that £ > 0, consistent with the fact that WIMPs 
close to the star are bound to it. 

For WIMPs on orbits that do not cross the star, the 
radial action equation can be solved analytically since 



C. Which WIMPs cross the star and when? 

For a WIMP to be captured, it has first to cross the 
star. Thus we need to follow how the WIMPs behave in 
and around the dark star. 

For a WIMP to be captured, it has first to cross 
the star. In the central regions of the halo, the star 
completely dominates the gravitational potential, as dis- 
cussed in the previous Subsection. The potential is thus 
given by '^^,{r), which is Keplerian outside the star and 
non-Keplerian inside. 

Once the star has formed, WIMP orbits not crossing 
the star are Keplerian. For these orbits, the minimal 
distance from the center of the star, rmin, can be written 
in closed form as a function of the adiabatic invariants 
and the mass of the star Af*. Using Eq. pO)) to find the 
smallest zero of the square root in Eq. ([7]) gives for this 
Keplerian case 



[Jr + 2njy 



(Jr + 2TTjy 



(11) 



This radius depends on time through the time depen- 
dence of the star mass Af^,. 

If rj^jj^ is greater than the radius of the star the 
orbit will not cross the star at this stage in the stellar 
evolution. Similarly, if r^^^ is smaller than the radius of 
the star i?^, the orbit will cross the star at this stage in 
the stellar evolution. One can be sure of these conclusions 
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FIG. 2: The evolution of Mt,Rt: during the star formation. 
First Mi,Ri, grows as the star grows but then drops fast as 
the star approaches the main sequence. 



because the Keplerian solution is unique, given E , J and 
the point mass potential. Hence, if the Keplerian orbit 
solution is found to be partly inside the star, no Keplerian 
solution with the orbit completely outside the star exists 
and so the true orbit must have parts inside the star. 
Conversely if the Keplerian orbit solution is found to be 
completely outside the star one has found the true orbit 
since the potential is Keplerian outside the star. 



As we are only interested in WIMPs on orbits cross- 
ing the star only orbits fulfilling rj^j^ < i?^, at a given 
time are interesting. As long as the individual WIMPs 
do not scatter, the adiabatic invariants in Eq. (fTT|) are 
conserved, implying that M^rJ^jj^ is constant in time (it 
is an adiabatic invariant). Along with the requirement 
that WIMPs crossing the star have to fulfill the relation 
-^*''min — Mi^Ri,, one finds that the star can potentially 
access the maximum amount of dark matter when Mi,Ri, 
is maximized. Fig. [5] shows the behavior of the quantity 
M-kRi, as a function of time for the Spolyar et al canoni- 
cal case. At first Mi^Ri, increases as the dark star grows, 
meaning that more and more WIMPs can cross the star. 
At the time of the Kelvin-Helmholtz (KH) contraction 
(t ^ 0.3 Myr), the quantity Mi^R^, decreases rapidly, leav- 
ing behind WIMPs that were once in stellar contact. For 
example, WIMPs with M^r|^^i„ - 10^ MqRq cross the 
star for times between 0.02 Myr and 0.3 Myr, but can 
no longer cross the star or be captured after 0.3 Myr. 
This makes it hard for the star to capture WIMPs in the 
later stages after the KH contraction. As a consequence, 
the accessible amount of dark matter that has survived 
the formation of the star might not be high enough to 
support the dark star after its KH contraction. 



D. The fate of the WIMPs 

Once we know which Monte-Carlo WIMPs cross the 
star at a given stage of the dark star evolution, the next 
step is to find their orbits. Some of these orbits are 
only partially inside the star, some are completely in- 
side the star. Knowing the gravitational potential inside 
and outside the star and the adiabatic invariants and 
J7r, enough information is given to determine E from Eq. 
([7]). This is done as explained at the end of Section II. B. 
Having determined the energy and the angular momen- 
tum of the WIMP in orbit, the full orbit is specified up 
to the spherical symmetry (that is up to the specification 
of the plane of the orbit and of the angle of the periilun 
[23). 

In the Monte-Carlo, to find the fate of the WIMP 
we need to calculate the probability per unit time for 
the WIMP to scatter off nucleons in the star as well as 
the probability for the WIMP to annihilate with other 
WIMPs in the star. In the following, we derive the scat- 
tering and annihilation rates at a given stage of the stellar 
evolution. 

When on an orbit crossing the star, the WIMP spends 
some time 



dr 



(12) 



in the shell between r and r + dr. Here the factor of 
2 accounts for the fact that the WIMP can cross the 
shell twice, once on its way in and once on its way out; 
moreover. 



Vr = - {Jlr? (13) 
is the (magnitude of the) radial velocity, where 



(14) 



is the speed at position r. The probability that in one 
passage through the star the WIMP annihilates against 
another WIMP in the dark matter population within the 
star is then given by 



^ann 2 



0'annt"^DM(?') , (15) 



where for the number density of WIMPs nDM(?') at ra- 
dius r we take the results from the dark star evolution 
in Ref. [2]. For WIMPs on orbits completely inside the 
star, the phrase "one passage though the star" refers to 
one revolution of the WIMP orbit. The number of times 
a WIMP passes through the star in the time At is 



At 



(16) 



where the orbital time T is determined numerically from 
the integral 



dr 

Vr 



(17) 
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Finally the probability that a WIMP annihilates during 
the time interval At is 

-Pann rji ^ann ■ (-^^) 

We annihilate Monte-Carlo WIMPs according to this 
probability. 

When passing the stellar material there is also some 
probability for the WIMP to scatter off nuclei in the star. 
Here only scattering off hydrogen is considered since we 
assume the scattering cross section is spin- independent, 
and so no scattering occurs off helium. Heavier elements 
are a negligible component of the material forming the 
first stars. We assume that at all radii hydrogen con- 
tributes 75% of the stellar mass density. 

In one passage through the star the probability for the 
WIMP to scatter in a shell at radius r is given by 

dr 

dT'scatt = 2 CTscatt V nuir) — . (19) 

Vr 

Here the factor of 2 accounts for the the fact that the 
WIMP crosses the shell twice, and 71h('") is the number 
density of hydrogen atoms in the dark star (taken from 
Ref. |3]). The scattering probability per WIMP passage 
through the star is then given by 

/'min(rmax ) 

^scatt 

/ dPscatt. (20) 

Then the probability that a WIMP scatters in the time 
At is 

_At 

-^scatt — rp "^scatt • v^-i) 

In a scattering event, the WIMP loses an amont of 
energy A£ subject to the condition 

A8<(!^^2^ymr)-S]. (22) 

Here \I'(r) — £ is the WIMP kinetic energy before the 
scattering at the scatterer location, and mn is the mass 
of the target nucleus, which in this paper is always hydro- 
gen since we assume the spin dependent scattering cross 
section to dominate. All possible values of the energy loss 
AE are equally probable, making it straightforward for 
us to determine the WIMP energy £' after the scattering. 

At a given evolutionary stage of the star, we use the 
scattering probability per unit time in Eq. (|2ip to decide 
how long it takes before the WIMP scatters. Then we 
compare this time to the time left in the given evolu- 
tionary stage to decide if the WIMP scatters before the 
dark star evolves to a new stage. Afterwards, the Monte- 
Carlo uses the annihilation probability per unit time in 
Eq. p8|) to evaluate if the WIMP annihilates before it 
scatters or before the current evolutionary stage ends. If 
the WIMP scatters, the Monte- Carlo uses the distribu- 
tion dT^scatt /fscatt to determine in which shell the WIMP 



scatters. In the scattering, the WIMP loses energy and 
ends up on a lower energy orbit, specified by new values 
of £' and J''. The new direction of the scattered WIMP, 
and hence J'', are easily determined as the scattering is 
isotropic in the center of mass system. After the WIMP 
is scattered, it is allowed to scatter again in the time left 
before the next evolutionary stage of the dark star. 

When passing to the next evolutionary stage, the adia- 
batic invariant JT/ is needed and can be determined using 
Eq. ([7]) with the £' and of the last orbit of the scat- 
tered WIMP. This allows us to follow the WIMP response 
to the stellar evolution whether the WIMP scatters or 
not. If the WIMP scatters many times it eventually sinks 
towards the core of the star and annihilates with the cen- 
tral WIMP population, injecting its energy into the dark 
star. 

Ref. [s^ ] contains further details and discussions on 
the scattering process. There one can also find how to 
determine the angular momentum of the WIMP after it 
has scattered. 



III. RESULTS 

Our Monte-Carlo follows a high number of WIMPs 
from the initial NEW halo until they annihilate or the 
star has formed. It allows us to track how much energy 
dark matter deposits into the star at the different stages 
of stellar evolution. Eig. [3] shows the resulting annihila- 
tion luminosities with and without scattering. The top 
solid (red) line includes WIMP-nucleon scattering, while 
the bottom dashed (blue) line does not. At the early 
stages of star formation, scattering of WIMPs off star 
nuclei is not important because the baryon density inside 
the dark star is low. As the dark star becomes denser, 
scattering becomes important, provided the scattering 
cross section is large enough. Eor the stellar evolution 
in [5| , the dark star becomes dense enough for scattering 
to be important just as it starts its rapid KH contraction 
onto the main sequence (see Eig. ^ . With the rapidly 
increasing stellar density at this stag,e capture becomes 
very efficient and via repeated scattering the captured 
WIMPs rapidly sink toward the center of the star where 
the dark matter density is high. This makes the captured 
WIMPs annihilate fast and gives the burst in luminosity 
seen in Eig. [31 As long as capture is efficient, the high 
central WIMP density is rapidly refilled with more cap- 
tured WIMPs. 

However, the burst lasts a short time because (1) the 
KH contraction is rapid and the volume of dark mat- 
ter available for annihilation decreases quickly, and (2) 
it takes a short time to use up all the WIMPs that can 
be efficiently captured by the dark star, and no other 
WIMPs are available for efficient capture. The latter rea- 
son ("death by starvation") is the most important one, as 
discussed in connection with Eig. [HI in Section IIVI The 
star is very fast at eating the WIMPs that it can eat eas- 
ily. After the contraction, WIMP annihilation continues 
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FIG. 3: The time evolution of the dark star luminosity from 
WIMP-WIMP annihilations. The top solid (red) line includes 
WIMP-nucleon scattering in addition to annihilation; the bot- 
tom dashed (blue) line does not include scattering. 
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FIG. 4: The dark matter mass accessible to the star as a 
function of time. In other words, the total mass of the WIMPs 
on orbits intersecting the star at the different stellar evolution 
stages. The dotted (blue) line includes both WIMP scattering 
and annihilation. The dashed (green) line includes WIMP 
annihilation but no scattering. The solid (red) line includes 
neither scattering nor annihilation. 



but at a much smaller rate; in the case we study, the anni- 
hilation after contraction is suppressed by (very roughly) 
about 15 orders of magnitude (five orders of magnitude 
coming from the change in stellar volume and two times 
five orders of magnitude coming from the decrease in den- 
sity due to annihilation and scattering, extracted from 
Fig. [5] below) . 

When WIMPs scatter, they lose energy, go on smaller 
orbits, and spend more time near the center of the star 
where the density is higher and the annihilation faster. 
Both annihilation and scattering remove WIMPs from 
orbits that cross the star. This is illustrated in Fig. 31 
where the mass of dark matter accessible to the star (i.e. 



FIG. 5: The density profile of dark matter around the star 
just before (dotted blue line) and just after (solid red line) 
the Kelvin-Helmholtz contraction at the end of the dark star 
phase (scattering is included when generating this graph). 
The vertical straight line marks the radius of the star af- 
ter contraction; before contraction, the radius of the star is 
4.2x10^^^ cm, about four times the length of the horizontal 
axis. 



the total mass of WIMPs on orbits crossing the star) is 
plotted as a function of time. The three lines illustrate 
the effect of the stellar evolution as well as that of scat- 
tering and annihilation. The top (solid red) line is similar 
to the curve in Fig. [51 and includes just the gravitational 
response of the WIMPs to the star and neither scatter- 
ing nor annihilation. It shows that the mere contrac- 
tion of the star already makes it a much smaller target 
for WIMP capture, by about one order of magnitude in 
the accessible mass. Annihilation reduces the number of 
WIMPs that cross the star (green dotted line), as anni- 
hilation removes WIMPs from the system. The inclusion 
of scattering reduces this number further, which might 
at first not seem obvious. The reason is that WIMPs on 
orbits not crossing the star cannot scatter in the star. 
Hence the total mass of WIMPs on orbits crossing the 
star can never increase due to scattering. Furthermore, 
scattering puts WIMPs onto lower energy orbits, making 
them more vulnerable to annihilation. Thus the num- 
ber of WIMPs in stellar contact is further reduced (see 
Fig. 13]). As discussed in Section lTVl our Monte-Carlo only 
includes WIMPs from the inner 1 percent of the NFW 
halo scale radius rs , but WIMPs from regions further out 
do not contribute significantly to the annihilation rate in- 
side the dark star. 

From the Monte-Carlo one can also extract the dark 
matter density profile around the formed star. This is 
shown in Fig. O for two different times, just before and 
just after the KH contraction (blue dotted line and red 
solid line, respectively). The solid red line is for the time 
0.4 Myr, which is just after nuclear burning has started. 
The blue dotted line shows the time 0.28 Myr, which is 
just before the star contracts enough for scattering to 
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be important, as can be seen in Fig. [31 The vertical, 
black line marks the size of the star at t = 0.4 Myr. At 
t = 0.28 Myr the stellar radius is 4.2 x 10^'^ cm and hence 
outside the graph. The unevenness of the graph comes 
from having simulated a finite number of WIMPs; the 
very central region has far less simulated WIMPs (as the 
fraction of WIMPs reaching so far in is very low) and 
hence greater inaccuracy. 

The KH contraction of the dark star increases the dark 
matter density in the central region through adiabatic 
contraction. The density peak at ~ 2 x 10^^ cm on the 
solid line in Fig. [5] is a residual of this density increase. 
Moreover, scattering increases the density at the very 
center by moving WIMPs into low energy orbits on which 
they spend more time near the center. The underdense 
"hole" at small radii on the solid line in Fig. [5] is a conse- 
quence of the rapid WIMP scattering and annihilation. 
WIMPs with large contributions to the dark matter den- 
sity inside the dark star are very vulnerable to scattering 
which makes them sink to the core and annihilate away 
fast. In the absence of inhomogeneities, this final config- 
uration is stationary: WIMPs in orbits that crossed the 
star have scattered and annihilated away, and WIMPs 
in orbits that do not cross the star remain in those or- 
bits. Even if the star would move in space, the WIMPs 
surrounding it would move with it, much like the Moon 
moves with the Earth while orbiting it. As long as this 
setting is not disturbed by inhomogeneities, the WIMP 
density profile around the star will look essentially the 
same until the star dies (except that the annihilation 
timescale near the density peak is short and annihila- 
tion will reduce the density). The star will then simply 
be a normal massive star, living its life undisturbed by 
its envelope of dark matter. 



IV. DISCUSSION 

We mentioned in Section IIIII that the death of a dark 
star is due to dark matter exhaustion rather than rapid 
KH contraction. In the first case, the WIMP supply runs 
out first and the lack of luminosity from WIMP annihi- 
lation causes the dark star to contract onto the main 
sequence. In the second case, the dark star contracts 
first, for example because the WIMPs do not provide 
enough luminosity, and is then unable to continue cap- 
turing enough WIMPs to remain a dark star because 
it has become too small. There is a third possibility, 
namely that the dark star does not contract or might 
even grow, since the extra luminosity provided by the 
captured WIMPs make it puffier and it may become more 
massive by accreting gas for a longer time. If the dark 
star does not grow, there is a maximum amount of WIMP 
mass that can be captured before going onto the main 
sequence, given by the graph in Fig. 2] at the time just 
before contraction ~ 0.27 Myr. 

Here we argue that even if the dark star does not con- 
tract, the supply of dark matter is rapidly exhausted (the 
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FIG. 6: The energy injected into a dark star from WIMP an- 
nihilation for a star with modified stellar evolution so that it 
is artificially frozen just before contraction. The WIMP lu- 
minosity in the altered stellar evolution is given by the solid 
(red) line; for reference, the WIMP luminosity without alter- 
ing the stellar evolution is shown by the dashed (blue) line. 
Both luminosities include scattering. This is to illustrate that 
the disappearance of dark matter in contact with the star is 
due to the depletion of star-crossing orbits with low angular 
momentum and is not an artifact of analyzing a star that 
contracts too rapidly. 



first case above). In principle, we should recompute the 
evolution of the dark star in the light of the capture re- 
sults we find, instead of using the dark star evolution 
tracks without capture in f?! (notice that we cannot use 
the evolution tracks with capture in [2] because they are 
based on Gould's capture formula that does not apply to 
our case). A detailed numerical study of the response of 
the dark star to the extra luminosity produced by cap- 
tured WIMPs is under way |28|. Here we avoid going 
into the details of the stellar response; we instead artifi- 
cially freeze the evolution of the dark star on the verge 
of contraction. We prolong the stage in the stellar evo- 
lution when the dark matter capture via scattering is at 
a maximum (the peak in Fig. [3]) from 6,000 to 60,000 yr. 
The dark matter luminosity for the star with this slight 
change in the stellar evolution is shown in Fig. |6l which 
is analogous to Fig. [3] but for the KH contraction which 
is artificially delayed as just explained. We still observe 
a drastic reduction of the luminosity due to WIMP anni- 
hilation. We find that the extra dark matter supplied via 
WIMP nucleon scattering can only support the star for a 
limited period of time; it does not result in a long lasting 
dark star phase. In other words, we conclude that the 
demise of the dark star phase is due to the exhaustion of 
the dark matter supply. 

In our code we have not included WIMP annihilations 
that occur outside the star. We expect these annihila- 
tions are not important for the energy production in the 
dark star for the following reason. The dark star can- 
not absorb the energy released by WIMP annihilations 
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outside of it, but the outside WIMP annihilations would 
somewhat reduce the number of WIMPs available for an- 
nihilation inside the dark star. However, during the dark 
star phase only a very small fraction of the accessible 
WIMP population would have annihilated away outside 
the dark star, making this depleting effect unimportant 
for the dark matter heating of the dark star. Dark mat- 
ter annihilations can also affect the cooling of the dark 
star by ionizing the gas (discussed in |2l|). We are here 
only looking at the dark matter component of the dark 
star process so feedback through gas ionization is beyond 
the scope of our paper. 

The conclusion that the main sequence star does not 
significantly benefit from its surrounding dark matter 
halo is based on the assumption that the final dark mat- 
ter configuration, as shown in Fig. [SJ is stable enough 
against perturbations to largely sustain its shape over 
time. For example a close encounter with another star 
could severely affect the dark matter envelope. How- 
ever, without long term heating from WIMP annihila- 
tions, stars this massive are very short lived, with ex- 
pected lifetimes of less than 2 Myr [33]. Since we assume 
there is only one dark star per dark halo, the time scale 
for the dark star to encounter another star is expected to 
be much longer than the star's lifetime. 

The star locally completely dominates the gravita- 
tional potential so if the star experiences a small acceler- 
ation towards some distant object the dark matter enve- 
lope would just follow, keeping its shape and distance to 
the star. If the star experiences a more violent event, so 
that it is suddenly kicked out of its position and thrown 
into its dark matter envelope, the star will once again 
be able to benefit from the dark matter. However, the 
system will fast reconfigure itself around the new grav- 
itational center at the stellar location. The star would 
quickly eat up the WIMPs it can eat easily, and the other 
WIMPs will have too high an angular momentum rela- 
tive to the new gravitational center to be reachable by 
the star. Such a drastic event would hence induce a short 
dark matter burst similar to the capture burst in Fig. [31 
followed once again by starvation of the star. For the star 
to have a long-term benefit from the dark matter enve- 
lope, continuous perturbations of the star or the dark 
matter envelope are needed. 

The initial star-forming gas cloud could also have some 
initial velocity relative to the dark matter halo. The 
forming star will still create a dark matter envelope 
around it. The star dominates the gravitational poten- 
tial and carries the WIMPs with it, hence the WIMPs 
in the star's dark matter envelope do not feel that the 
halo has a net velocity relative to the star. In brief, the 
initial conditions of the WIMP orbits in the dark matter 
envelope will be slightly different but the dynamics of the 
system will be the same, so the final result for the dark 
matter in the central region of the halo will essentially be 
the same. The only WIMPs that will notice the relative 
velocity between the star and the halo are the WIMPs 
traveling far enough out in the halo to also probe the 



halo's gravitational potential. However, as these WIMPs 
pass the region where the star is, their velocities are very 
high, making any reasonable stellar velocity negligible in 
comparison. We also show below that these WIMPs are 
not important anyway. Finally, a net velocity of the star- 
forming cloud could have a significant impact only if the 
star would form far away from the center of the halo and 
thus miss the halo central high dark matter density. As 
long as the formation is not very far from the center, the 
off-center location would primarily reduce the initial dark 
matter heating from adiabatic contraction as this con- 
traction starts before the star has taken its dominance 
over the local gravitational potential. 

The star's host dark matter halo is expected to have 
substructure with smaller subhalos of dark matter. The 
effect of these subhalos on the formation of dark stars 
is still to be explored. In any case, subhalos in the cen- 
tral region are likely to be destroyed by the formation 
of the dark star. Since subhalos are only a fraction of 
the total dark matter population around the star, and 
since each subhalo would provide WIMPs for a limited 
amount of time and there would be gaps in the WIMP 
supply between subhalo crossings, the subhalos are not 
expected to sustain a dark star through capture over a 
long period of time. It could be possible, however, that 
large- and small-angle gravitational scattering of WIMPs 
off subhalos around the dark star changes the phase-space 
structure of the dark matter envelope. A similar effect 
would arise from inhomogeneities in the gas surrounding 
the dark star. Many small-angle deflections could over 
time make some WIMPs in the envelope end up crossing 
the star. 

The part of the dark matter envelope closest to the 
star should be the most vulnerable to perturbations that 
make the WIMPs cross the star. If the perturbations 
would be such that the star first accesses the WIMPs in 
its immediate surroundings, it would need to capture all 
of the dark matter out to a radius of lO^i?* (after KH 
contraction) to prolong its lifetime by roughly 1 Myr. 
Hence, even if the dark star could efficiently capture 
the dark matter in this closest reservoir, it is still not 
enough to substantially prolong the life of the star. To 
do so, the dark star needs to efficiently access dark matter 
from outside lO^i?^. By a simple geometrical argument, 
WIMPs further away from the star should be less likely 
to be perturbed by inhomogeneities into orbits crossing 
the star. Orbits stretching further out also have longer 
orbital times, giving fewer stellar crossings per unit time. 
Assuming that all WIMPs are constantly subject to per- 
turbations, and using a simple geometrical probability 
factor [Ri^/rY for the probability per orbital revolution 
that an orbit extending out to radius r intersects the star, 
gives an estimate of the capture rate and of the luminos- 
ity contribution from WIMPs on orbits stretching further 
out than lO'*/?^. This luminosity is too small to sustain 
the dark star phase. 

The overall conclusion of this is that it could in prin- 
ciple be possible for WIMPs from the dark star's WIMP 
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envelope to temporarily feed the main sequence star with 
enough dark matter to relaunch a short dark star phase, 
for example through perturbations of or direct interac- 
tion with dark matter subhalos in the dark matter halo. 
What the above argument shows is however that it is very 
unlikely for the dark matter envelope to be able to con- 
tinuously supply a dark matter capture rate high enough 
to obtain a dark star phase comparable or longer than 
the star's main sequence lifetime. 

In the tradition of previous papers about Pop III dark 
stars, this work assumes that the first stars form individ- 
ually near the center of a dark matter halo. Some recent 
work [s^l has suggested that some of the first stars might 
have formed as binary systems. If that would be the case, 
the analysis of the dark matter response during their for- 
mation would be severely more complicated. In general, 
one would expect a less important density increase due 
to adiabatic contraction and also, once the nearby star 
has formed, their motion should accelerate and throw out 
WIMPs bound to the central region, further reducing the 
central density. To make this statements more precise, a 
specific investigation is needed. 

Population III stars as massive as discussed here can 
become unstable to pulsations [ssl, fssj. An interesting 
thought is that as the star pulsates it has access to a part 
of its dark matter envelope. This could potentially have 
interesting consequences, but would require a study in 
itself. In any case, the dark star can only eat up the dark 
matter in the envelope during the first few pulsations, 
since in subsequent pulsations the dark matter supply 
will have been depleted. Moreover, the gas density in 
the outer layers of the pulsating star may be too small 
for scattering of WIMPs to occur with any significance. 

To speed up the Monte-Carlo calculation, when pick- 
ing WIMPs from the initial NFW halo, we restrict our- 
selves to the region r < COlr^. In principle, one should 
also consider WIMPs further out, but the probability 
that they scatter in the dark star decreases rapidly with 
distance and they soon become irrelevant. In addition, 
the orbital time for WIMPs with orbits extending out to 
0.02rs is ~ 1 Myr, meaning that these WIMPs will not 
be captured by the dark star within the much shorter 
lifetimes we find. We have checked that WIMPs beyond 
O.Olrg are indeed irrelevant for our results by running 
a Monte-Carlo simulation picking WIMPs from a region 
that extended tenfold in radius, r < O.Irs. 

In our derivation, we have neglected the contribution 
from WIMPs that are not bound to the star. Here we 
present several arguments leading to the conclusion that 
the maximum rate at which WIMPs can be captured by 
a dark star from the unbound WIMP population is in- 
deed small enough to be completely negligible. Unbound 
WIMPs are on orbits that extend beyond the radius of in- 
fluence of the star rjnfl. This radius was estimated in Sec- 
tion II. B to be Tinfl ~ O.OIts for the capturing stages of a 
dark star, namely the final stages with Mi, ^ 10^ Mq and 
i?* = 4x 10^^ cm. Now, the density of WIMPs at r rjnfi 
is comparable to the original density in the NFW halo, 



namely ^ lO^GeV/cm^. This density is much smaller 
than the density of bound WIMPs near the star, which 
is of order 10^^ GeV/cm'^. Therefore unbound WIMPs do 
not contribute significantly to the WIMP density near the 
star. One might worry that gravitational focusing would 
be able to collect more unbound WIMPs than the previ- 
ous argument would indicate. However, the focusing ef- 
fect can be included in an estimate of the amount of dark 
matter mass that could fall in from a distance of nnfi and 
cross the star. One finds for this rate of mass collection 

^unbound ^ T^RlF^i whcrC = F(rinfl)Wesc/'f^(nnfl)^ 

is the inward flux of unbound WIMPs at the surface 
of the star, Wcsc/^(nnfi)^ is the gravitational focusing 
factor, and F^ri^fi) = p(rinfl)i;(rinfl) is the inward flux 
at r = finfl. Estimating the density and speed of 
WIMPs ai r — rina using pnfw and f (rinfl) = 

[GA/NFw(^infl)/''infi]^^^ 1 km/s, givcs a mass collec- 
tion rate of ^ 10~^ Mo/Myr from the population of un- 
bound WIMPs. To survive at a luminosity of ^ 10^ L©, 
the star would need an inflow of dark matter WIMPs 
to annihilate equal to ~ IMQ/Myr, much higher than 
what unbound WIMPs could provide. Also, only a small 
fraction of the WIMPs passing through the star will ac- 
tually scatter and be captured by the star, reducing the 
heating from unbound WIMPs further. In addition, the 
orbital periods of WIMPs on unbound orbits, i.e. orbits 
that reach out to r = Tina or more, is at least 2 Myr (as- 
suming Keplerian orbits in this estimate). This time is 
comparable to the lifetime of the dark star, including the 
later main sequence stage. Thus unbound WIMPs can- 
not reach the star in time to keep it alive through dark 
matter annihilation. 

It has recently been argued in |[20,] that triaxiality of 
the dark matter halo could efficiently continuously sup- 
ply WIMPs to the star and maintain the dark star phase 
on cosmological time scales. However, as the star dom- 
inates the gravitational potential in the central part of 
the halo, only WIMPs traveling out to radii comparable 
to the star's radius of influence rjnfl feel the triaxiality 
of the dark matter halo. Thus the method in the previ- 
ous paragraph can be used to estimate the contribution 
to the capture rate from WIMPs probing the triaxiality 
of the halo. This estimate gives that the supermassive 
dark star's capture rate from the triaxial halo in [20] is 
too low by a factor of ^ 10^ to sustain the luminosity 
of the star. Even though this estimate is based on the 
outer halo being isotropic, accounting for anisotropics or 
chaotic orbits is not expected to increase the capture rate 
by a factor of ~ 10^. A proper treatment of the triax- 
iality may slightly increase the capture rate thanks to 
elongated orbits that are bound to the star and almost 
reach rinfl, but again the increase cannot be by 9 orders 
of magnitude. In summary, a triaxial halo could feed the 
star with dark matter for a very long time but as the cap- 
ture rate is too low to sustain the dark star, triaxiality 
of the halo is not interesting for the capture rate. 

An intriguing possibility is that long-lived dark stars 
may form from the ashes of Pop HI dark stars (Phoenix 
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Dark Stars). When the main sequence star following 
the dark star phase explodes, it ejects gas into a re- 
gion that can still have a considerable density of dark 
matter WIMPs. This dark matter is the remnant of the 
dark matter envelope that surrounded the dark star. In 
fact, for (TannV = 3 X 10^^^ cm'^/s, a remnant dark matter 
density uq — l/((Tann^'^o) ^ 10^ WIMPs/cm'^ can be sus- 
tained for a time Iq ^ 10 Gyr (approximately the age of 
the universe). For a 100 GeV WIMP, this corresponds to 
a mass density po — 10^°GeV/cm'^ = 2 x 10^^^ g/cuv^. 
Some of the ejected gas may fragment and form Pop II 
and later Pop I stars, which can move relative to the 
dark matter remnant and capture WIMPs from the envi- 
ronment in the way described by previous authors j^-d] ■ 
These "daughter" stars are expected to be less massive 
than their parent Pop III dark stars, and thus would need 
less dark matter to be a dark star. Given appropriate en- 
vironmental conditions, some of these second generation 
dark stars might perhaps survive until today. 

V. CONCLUSIONS 

As found by previous authors, the density of dark mat- 
ter inside a Pop III protostar can be high enough to lead 
to a dark star phase powered by dark matter annihilation 
instead of nuclear fusion. A large WIMP-nucleon scatter- 
ing cross section may further increase the energy injected 
into the dark star by WIMP annihilations. As discussed 
in the introduction, many authors have claimed that the 
WIMP supply from scattering may last a long time, even 
until the present day. Our work does not support those 
findings. 

We find that WIMP capture is not important until the 
dark star contracts on its way to the main sequence and 
becomes dense enough to be a viable dark matter target. 



We also find that WIMP capture lasts a relatively short 
amount of time as the WIMPs available for capture are 
quickly depleted in less than lO'^ years and their orbits 
are not replenished. 

It is also very unlikely that perturbations in the star's 
dark matter envelope could keep the inflow of dark mat- 
ter high enough to support a further prolongation of the 
dark star phase, namely for lifetimes much longer than 
the 10^ years we here calculated. Furthermore the inflow 
of dark matter from the outer regions of the halo is too 
small to be of any interest. From this we conclude that 
a single Pop III star forming at the center of an NFW 
halo will have a total lifetime of order a few million years 
and that the dark star phase supported by WIMP cap- 
ture will not substantially prolong the lifetime of the dark 
star. 

Our work indicates that WIMP capture in Pop III dark 
stars is not able to make them live forever. However, 
later generations of less massive stars produced in the 
dark matter environment left behind by Pop III dark 
stars may capture enough dark matter to form a second 
generation of dark stars. 
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